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We describe a method to express the susceptibility and higher derivatives of the free energy 
in terms of the scaling variables (Wegner's nonlinear scaling fields) associated with the high- 
temperature (HT) fixed point of Dyson hierarchical model in arbitrary dimensions. We give a 
closed form solution of the linearized problem. We check that up to order 7 in the HT expansion, all 
the poles ("small denominators") that would naively appear in some positive dimension are canceled 
by zeroes ("small numerators"). The requirement of continuity in the dimension can be used to lift 
ambiguities which appear in calculations at fixed dimension. We show that the existence of a HT 
phase in the infinite volume limit for a continuous set of values of the dimension, requires that this 
mechanism works to all orders. On the other hand, most poles at negative values of the dimensional 
parameter (where the free energy density is not well-defined, but RG flows can be studied) persist 
and reflect the fact that for special negative values of the dimension, finite-size corrections become 
leading terms. We show that the inverse problem is also free of small denominator problems and 
that the initial values of the scaling variables can be expressed in terms of the infinite volume limit 
of the susceptibility and higher derivatives of the free energy. We discuss the existence of an infinite 
number of conserved quantities (RG invariants) and their relevance for the calculation of universal 
ratios of critical amplitudes. 

PACS numbers: 05.10.-a, 05.50.+q, ll.10.-z, 64.60.-i 
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I. INTRODUCTION 

In many problems, one faces the challenge of deriving 
the macroscopic consequences of a microscopic theory. 
As we look at the problem at increasingly large scales, 
a sequence of effective theories appear and under some 
appropriate conditions, an infinite volume limit can be 
taken. A general method that allows us to construct 
these flows in the space of theories is the Renormaliza- 
tion Group (RG) method 0. The study of some RG 
fixed points and of the linearized flows close to these fixed 
points has produced a successful picture of the universal 
behavior in second order phase transitions. On the other 
hand, controlling the RG flows beyond the linearized ap- 
proximation and calculating the related nonuniversal be- 
havior are more difficult issues. This is unfortunately 
necessary to calculate the critical amplitudes. 

As a first step, one can deal with the nonlinear RG 
flows for simplified models where the RG transforma- 
tion can be implemented without major technical diffi- 
culty. One possibility is to use approximate versions of 
the exact RG equations 0, El such as the local potential 
approximation J5L Another possibility to address nonlin- 
earquestions [UjS0 is to use Dyson's hierarchical model 
In the following, we use this lattice model for which 
the block-spin method can be easily implemented. This 
model is briefly reviewed in section^ Other approaches 
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of nonlinear aspects of the RG flows can be found, for 
instance, in Refs. HH EJ H 03 

In the context of ordinary differential equations, a stan- 
dard method |l4j to go beyond the linearized approxi- 
mation in the vicinity of a fixed point, consists in con- 
structing a new system of coordinates where the equa- 
tions become linear. However, this type of procedure 
is often plagued with the "small denominator problem" 
initially encountered by Poincare in his study of per- 
turbed integrable Hamiltonians. In the context of the 
RG method, these new coordinates are called the scal- 
ing variables (or the nonlinear scaling fields) and were 
first introduced by Wegner 15]. Recently, we have pro- 
posed an ab-initio calculation of the critical amplitudes 
in the high-temperature (HT) phase of this model • In 
this calculation, the critical amplitudes are RG-invariant 
made out of the nonlinear scaling variables associated 
with Wilson's nontrivial IR fixed point and the nonlin- 
ear scaling variables associated with the HT fixed point. 
In this approach, the two fixed points are in some ap- 
proximate sense dual to each others. The scaling 
variables associated with Wilson's fixed point have been 
extensively discussed, but much less is known about those 
associated with the HT fixed point. We emphasize that 
being able to use both kind of variables is quite conve- 
nient for the study of the RG flows in the intermediate 
region between the two fixed points. 

At first sight, the construction of the scaling vari- 
ables associated with the HT fixed point is impossible 
for D = 3 and more generally for rational values of D, 
because some of the denominators are exactly zero. How- 
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ever, a numerical study in D = 3 showed 17] that in all 
of the 36 zero denominators considered, a zero numerator 
miraculously appears. This strongly indicates the exis- 
tence of a general mechanism enabling us to overcome 
the small denominator problem. 

In this article, we show that such a mechanism ex- 
ists and is closely related to the existence of the infinite 
volume limit of the susceptibility and higher derivatives 
of the free energy. In addition, we address the issue that 
whenever zero numerator and denominator appear at the 
same time, the coefficients of the nonlinear expansion ap- 
pear to be undetermined (jj). We show that this inde- 
terminacy can be lifted by a procedure similar to the 
dimensional regularization [Tsj used for the evaluation 
of Feynman diagrams. It should however be emphasized 
that it is not used here to take care of a UV problem, 
since we will be working with a lattice model. We will 
consider the construction of the HT scaling field in arbi- 
trary dimensions. In our construction, the zero denomi- 
nators appear as poles at particular dimensions and one 
can study the mechanism of cancellation close to a pole 
but not exactly at the pole. 

For Dyson's hierarchical model, the dimension D ap- 
pears in a continuous parameter c = 2 1 ~ 2 / D introduced 
explicitly in section [H] The infinite volume limit is well 
defined for < c < 2, or in other words D > 0. The 
linear variables associated with the HT fixed point are 
introduced in section IIIII A closed form expression for 
the linear transformation which diagonalizes the linear 
RG transformation is given in arbitrary dimensions. The 
restriction to the first l max of these variables can be in- 
terpreted as a HT expansion. In section llVl we expand 
the linear variables in terms of the scaling variables . We 
show that up to order 7 in the HT expansion, the poles 
corresponding to zero denominators in positive dimen- 
sions (0 < c < 2) are exactly canceled by a zero at the 
numerator. The coefficients of the expansion are then un- 
ambiguously defined rational functions of c with no poles 
for < c < 2. Their poles appear only at negative values 
of the dimension where the statistical mechanics model 
does not have a well defined infinite volume limit. 

The linear variables are linear combinations of the av- 
erage values of the total field J2 X § x ■ ^ n sec ti° n E3 we 
use this fact to reexpress the connected parts of the av- 
erage values of the total field divided by the volume, or 
in other words, the susceptibility and the higher deriva- 
tives of the free energy density, in terms of the scaling 
variables. We show that up to order 7, the linear con- 
tribution is the only leading term in the infinite volume 
limit. In section IVT1 we explain why this should happen 
to all orders. In section lYlIl wc explain why it guarantees 
the cancellations discussed in section Hvl to all orders. 

Having showed that it is possible to construct a so- 
lution of the RG flows in the HT phase, we then need 
to calculate the initial values of the scaling variables in 
terms of the local measure (for instance, a Ising mea- 
sure or a Landau-Ginzburg measure) used to specify the 
statistical mechanics model. This amounts to inverting 



the previous expansions. In section [Villi we construct 
the scaling variables in terms of the linear variables and 
show that the coefficients are free of poles for < c < 2. 
We also show that, up to numerical constants, the ini- 
tial values of the scaling variables are the infinite volume 
limit of the susceptibility and higher order derivatives of 
the free energy density. This concludes our construction 
of a complete solution of the RG flows in the HT phase. 
To be precise, we have shown that various expansions 
can be constructed order by order without encountering 
any small denominator problems and that it is possible 
to study empirically the convergence of these series. In 
section llXl we show with an example how everything we 
have done can be used to calculate the HT expansion 
at finite volume. We also check explicitly that it yields 
results in agreement with calculations performed using 
independent methods [l^. In section Ixl we discuss the 
existence of an infinite number of conserved quantities 
and their relevance for the calculation of universal ratios 
of critical amplitudes. 

II. DYSON'S HIERARCHICAL MODEL 

In this section, we remind some basic facts about 
Dyson's hierarchical model that will be needed in the 
following. For more details, the reader may consult Refs. 
pfl l2lT | . We consider fields located at 2™ max sites labeled 
with n max indices x nmax , x%, each being or 1. We 
divide the 2™ max sites into two blocks, each containing 
2"max-i s ites. If x„ max — 0, the site is in the first box, if 
i„ mai = 1, the site is in the second box and so on. The 
non-local part of the energy reads 



-, "max 

*>* - ~\ E (|r E ( E h 




n—l x nrnax ,...,x n +i x n ,...,x± 

The partition function for a constant source J (or exter- 
nal magnetic field) reads 

Z(J) = II / d^ x W(<f> x )eM~l3H nt +Jj2 ( l ) y) ( 2 ) 

x ■> y 

We call W{4>x)d<f> x the local measure. The most common 
examples are the Ising measure W{4>) — <5(</> 2 — 1) or the 
Landau-Ginzburg measure W{4>) — exp(— A(f> 2 — Btf) 4 ). 
The RG transformation consists in integrating over the 
fields keeping their sum constant in increasingly large 
boxes. A fter each integration the fields are rescaled by a 
factor yj c/4 in order to keep the form of H n i identical, 
and the RG transformation generates a flow in the space 
of local measures. 

Note that for a constant configuration where all the 
fields take the same value <fi, the nonlocal part of the 
energy takes the value 

H nl {4>) = -2— {4>f\ E(|) n ■ (3) 
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In the infinite volume limit (n max — ► oo), the sum con- 
verges only for |c| < 2. Proofs of the existence of the 
thermodynamical limit for a Ising measure |22l |23| re- 
quire that the energy does not scale faster than the num- 
ber of sites. This means |c| < 2 for the model considered 
here. 

In the following, we will make a change of variables in 
order to get rid of /3 in front of H n i in Eq. and reab- 
sorb it in the local measure. Our main object of study will 
be the generating function (obtained by Fourier trans- 
forming the local measure) 



R n (k) = 1 + a n:1 k 2 + a n , 2 k 4 + . . . , (4) 



with 



a n ,i = (-0)'^|(f) ,n < ( E ■ ( 5 ) 

2 n sites 

The RG transformation can be summarized in terms of 
the recursion formula 



Rn+i(k) = C n+ i exp 



1 d 2 1 

29P 



^~ck 



-i 2 



(6) 



We fix the normalization constant C n so that i?„(0) = 1. 

It is important to remember that in the notation a n> i, 
the first index refers to the number of RG steps and the 
second to the powers of the total field. Sometimes, the 
number of RG steps n will be omitted, sometimes the 
vector index I will be replaced by boldface notations. 
We use the parametrization c = 2 1 ~ 2 / D such that a free 
massless field scales in the same way as in a usual D- 
dimensional theory. For reference, Dyson's parametriza- 
tion [U, was c = 2 2 ~ a . The logarithm of R generates 
the connected zero-momentum Green's functions at fi- 
nite volume. We emphasize that in the following, the 
temperature dependence has been absorbed in the ini- 
tial Ro(k). For instance, in the case of an Ising measure, 
R (k) = cos(V/3fc)- 

In the HT phase, polynomial truncations of order L 
in k 2 provide rapidly converging approximations 0, [2C 
The RG flows can be expressed in terms of a quadratic 
map in a l max dimensional space 



U n ,l 
Un,0 



with 



and 



rr = (c/4) 



(-i/2y+— (2(,i + v ))\ 



(7) 



(8) 



(9) 



for \l + v > a and zero otherwise. We use "relativistic" 
notations. The greek indices /i and v go from to l m ax, 
while latin indices i, j go from 1 to l max . Repeated in- 
dices mean summation unless specified differently. With 



the normalization of Eq. J7J), a„,0 = 1 for any n and 
is not a dynamical variable. Note that a truncation to 
order l ma x is always implicit in the following. However, 
for reasons that will be explained in the coming section, 
there is no explicit dependence in l m ax- 



III. THE LINEAR RG TRANSFORMATION 

In this section we discuss the linearized RG transfor- 
mation near the HT fixed point ai = for all i < 1. For 
small departure from the HT fixed point 5a n ,i the linear 
RG transformation reads 



Sa n +i,i — MlSa n , , 



with 



M * 4 V ( 2> (2i)\(j-i)\ 



(10) 



(11) 



for i < j and zero otherwise. 

The diagonalization of M. is not too difficult because 
of its upper triangular form. The spectrum is given by 
the diagonal elements: 

A (r) = 2(c/4r , (12) 

in agreement with Ref. @- We need to construct 1Z, a 
matrix of right eigenvectors, such that 

\r)l^l , (13) 



M\Tl\ = X (r) TZJ , 



(with no summation over r). For convenience, the 
columns of 1Z are ordered as the eigenvalues, < c < 4 
being assumed. We will then introduce the linear coor- 
dinates h ni i defined by 



and which transform as 



' / ri-\-1,r 



(14) 



(15) 



in the linear approximation. The matrix IZ^ and its in- 
verse are also upper triangular. This implies that h n ^ 
is of order f3 , just as a n ,i is. We will fix the normal- 
ization of the right eigenvectors in 1Z in such way that 
all the diagonal elements are 1. This guarantees that 
K,i =a lM + 0(/3 /+1 ). 

Before entering into the technical details of the con- 
struction of 7Z, an important consequence of the upper 
triangular form of M. should be noticed. The eigenvec- 
tors and eigenvalues of M. are independent of a possible 
truncation. In other words, the fact that 1Z is upper tri- 
angular means that the polynomial truncations of R to 
order k 2lma3: mentioned in section^ are indeed a projec- 
tion in the subspace spun by the first l max eigenvectors 
of M. 

We now construct 1Z. We first notice that for j > i, 



K 



o _ 



"4-c 



(16) 
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with V\ c-independent For indices no larger than 7, the 
entries of P are 



(I 6 


45 


420 


4725 


62370 


945945 \ 


1 


15 


210 


3150 


51975 


945945 





1 


28 


630 


13860 


315315 








1 


45 


1485 


45045 











1 


66 


3003 














1 


91 


\0 














1 / 



The matrix P has remarkable properties 



p? +1 p. 



i+l 
-l-ni+3 
' i+l 



= IP, 



i+2 



-pi+2-pi+3 _ <ni+l<r>i+2<r>i+3 
1 i 1 i+2 ~ 1 i 1 i+l 1 i+2 



(17a) 
(17b) 



and higher order ones that as we will see are related to the 
very simple form of the inverse matrix. We can condense 
these relations into the more compact recursion 



V 



i+mjyi+m+q 

i+m 



m + q 

q 



yi+m+q 



This implies the closed form expression 



1 { r (2i)!(j- 



0" 



(18) 



(19) 



Using Eq. H16|l and l|19|l . it is easy to check that, as a 
consequence of the binomial formula, we have provided 
an exact solution of Eq. i|13|) with the required normal- 
ization (ones on the diagonal). Similarly, one can show 
that the inverse has the very simple form 



(20) 



This equation does not hold for an arbitrary upper tri- 
angular matrix. It implies the identities l|17|) and many 
others. 



IV. EXPRESSION OF THE LINEAR 
VARIABLES IN TERMS OF THE SCALING 
VARIABLES 

In this section, we express the linear variables hi in 
terms of the (non-linear) scaling variables yi for which the 
approximate multiplicative transformation of Eq. H15|) 
becomes exact. If we use ln(yi) as our new coordinates, 
the RG flows become parallel straight lines. All the dy- 
namics is then contained in the mapping that we now 
proceed to construct. 

We first rewrite the RG transformation in the hi coor- 
dinates. Starting with the basic Eq. (JJJ, we replace ao 
by 1 and ai by nfh p . This yields 



h n +i,i — 



A(i)/l n ,i + h nt ph njq 



1 



h n p 



A pq h h 



(21) 



with coefficients calculable from Eq. ©. For instance, 



In general, upper roman indices transform with n and 
the lower ones with (7?.) -1 . By construction, the linear 
transformation is diagonal. 

We then introduce the expansion 



hi =yi+ 



s l,iii2- 



.VlV2 



(22) 



where the sums over the i's run from to infinity in each 
variable with at least two non-zero indices. In the follow- 
ing, we use the notation i for (ii, . . ). More generally, 
vectors will be represented by boldface characters. The 
unknown coefficients s/4 in Eq. 1221) are obtained by 
matching two expressions of h n +\ t i, one obtained from 
the RG transformation of the hi given in Eq. I|21() , the 
other obtained by evolving the scaling variables accord- 
ing to the exact multiplicative transformation 

Hn+1,1 = \l)Vn,l ■ (23) 

The matching conditions can be expressed as: 



/i„+i,i(h n (y)) = h n ,i(\iy n ,i, 
and yield the conditions 



A22M,: 



) 



(24) 



(25) 



with 



+ Ej+k+r=i s 'jn m A|™ ) Ag' 3 ,sp, k s 9 , r . (26) 

and 

A,i=A (() -nt) ■ (27) 

rn 

For a given set of indices i, we introduce the notation 



(28) 



One sees that Tq is the degree of the associated product 
of scaling variables and X\ its order in the HT expansion 
(since yi is also of order Given that all the indices are 
positive and that at least one index is not zero, one can 
see that if j + k = i then T q (j) < T q (i) and I 9 (k) < I q (i)- 
Consequently, Eq. I|26|l yields a solution order by order in 
Zo or in X\ (since the r.h.s is always contains s;^ of lower 
order in Tq or Ti ) provided that none of the denominators 
are exactly zero. The main goal of this article is to 
investigate what happens when some of the denominators 
happen to be exactly zero. 

Using the explicit expression of the eigenvalues Eq. 
(|12f) . we can rewrite the denominators as 



A,i = 2(- 



2 J (i)(-f:\Ii(i) 



(29) 
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Using the parametrization c = 2 1 2 / D , the zero denomi- 
nators appear when 

D -l(D + 2) = £>2o(i) -(£> + 2)Xi(i) . (30) 

Given that the X ? are integers, this can only occur at 
some rational values of D. Ignoring temporarily this set 
of values, we can say that for generic values of c, the 
denominator is not zero. Following the basic idea of di- 
mensional regularization, we will then perform, order by 
order in I\, the construction of the si^ for a generic value 
of c and discuss the limit where c takes some special value 
at the end of the calculation. 

We now determine the range of values of Tq and I\ 
relevant for our problem. In Eq. I|22l) . we have assumed 
that hi ~ yi for sufficiently small values of the scaling 
variables. The linear problem is completely solved and 
we may assume Xo(i) > 1. In addition, since both hi and 
yi are of order we needli(i) > /. At lowest non-trivial 
order in /3, we have = I, and 

D lti = {- A )\2-2 X °V) . 

In this special case, the only possible poles are at c = 0. 
However, the factor (|) ! at the denominator is exactly 
canceled by the same factor appearing in the A^ q in Eq. 
1)21 [I. More precisely 

h n+u = (|)'(2ft n ,, + ]T h n , p h n , q ) +0((3 l+1 ) . (31) 

p+q=l 

Using this, it is not difficult to prove by induction that if 
X(i) = Z, 

-m = n A (32) 

m 

It is thus clear that at the lowest non-trivial order, the 
coefficients have no singularities. 

We now discuss the case > I. We have in general 

Di, i = 2(^ l (c c „ t .) 1 - 1 ^ Ti, h (33) 

with 

Ti,i = (&r 1 -^- 1 ), (34) 

and 

c m , = 4x2( 1 - I »( i »/( I '( i »^ . (35) 

One should always keep in mind that c cr u. is a function 
of both I and i. Inspection of Eqs. 1© and 116(1 shows 
that the numerator has a factor c^Wfc — 4)'- I i( 1 ). Con- 
sequently, 

ai,i = (^~;) Il(i) - l Qu(c) , (36) 
c — 4 



where Qz,i(c) is a rational function of c with no poles or 
zeroes at or 4. We do not have a compact formula for 
these rational functions, however it is easy to calculate 
them using symbolic manipulation programs. 

Naively, we would expect that Qz,i(c) has a factor 
Ti t i(c) at the denominator and other poles inherited from 
the si y i of lower orders. The values of Qy(c) up to order 
/3 4 are shown in Table [I] The naive expectations con- 
cerning the poles are only observed in 9 cases out of the 
17 considered. In the 8 other cases, some cancellations 
occur. For instance, there is no (c+ 2) at the denomina- 
tor of Qi,(3,o,...)' More importantly, whenever c cr it. < 2, 
we observe a cancellation of all the factors appearing in 
Ti^c). This occurs, for instance, for Q2,(3,o,... )> where the 
factors c — 1 cancel. If we do the calculations explicitly 
using Eq. 1(26(1 . we obtain five terms at the numerator: 

— 11c 3 3 c 3 

^2,(3,0,...) = -^- ~ 8( _ 4 + c) 

c 3 15 c 4 

+ 4 (-4 + c) (-2 + c) + 64 (-4 + c) 



8 (-4 + c) (-2 + c) ' 
while the denominator reads: 

c 2 

A2,(3,o,...) = y( c ~ !) 

After reduction and factorization, the numerator be- 
comes: 

(-l + c)c 3 (6 + c) 

iV2,(3,0,...) lg ( _ 4 + c) (_ 2 + c) ' 

canceling the c — 1 at the denominator. We have pursued 
the same procedure up to order /3 7 and considered the 175 
possible terms. In 50 cases, we had c cr n. < 2. In each 
of these 50 cases, we observed a complete cancellation of 
T; ; i(c). It seems thus reasonable to conjecture that Qi t \{c) 
has no poles for |c| < 2. If this conjecture is correct, 
dimensional regularization provides a unique continuous 
expression for the coefficients for any c with |c| < 2 and 
the model is "solvable" using the recursion for the coef- 
ficients given by Eq. ((26(1 . Note that for values of c real 
and positive, the correspondence c = 2 1 ~ 2 / D implies that 
the interval < c < 2 corresponds to < D < +oo. The 
conjecture implies that for any value of c in this inter- 
val, we can construct analytical expression of a n j (which 
contains all the thermodynamical quantities) in terms of 
ao ; (which depends on the initial energy density): 

a n ,i = {K- l )thrMvi(aa), A£y 2 (a ), . . . ) . (37) 

The initial values of y(ao) have a simple interpretation 
discussed in section IVIIII 
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TABLE I: Values of Qi,i(c), c crit . and Tbfc) defined in the text. 



1 


Ilm 2/m 


Qu(c) 


Ccrit. 


-T M (c) 


1 


J/1 2 


z + c 
-2 + c 


2 


-2 + c 


1 


J/1 3 


— (4 — 20 c 4- c 2% l 
2 (-2 + c) 2 


2 


-4 + c 2 


1 


J/1 J/2 


-3 (-40 + c 2 ) 
-8 + c 2 


2^2 


-8 + c 2 


1 


J/1 4 


-120 + 156 c- 18 c 2 +c 3 
2 (-2 + c) 3 


2 


-8 + c 3 


1 


2/1 2 2/2 


3 (-11520 - 640 c + 1184 c 2 + 288 c 3 + 40 c 4 - 26 c 5 + 3 c 6 ) 


223 


-16 + c 3 


(-2 + c) (-8 + c 2 ) (-16 + c 3 ) 


1 




1536 
-32 + c 3 


2 

223 


-32 + c 3 


1 


m 2/3 


15 


223 


-32 + c 3 


2 


2/i 3 


6 + c 
2 (-2 + c) 


1 


-1 + c 


2 


2/i 2/2 


1 /I _L ^ 
14 -+" C 

-2 + c 


2 


-2 + c 


2 


4 

2/1 


-(-44- 28c + c 2 ) 
4 (-2 + c) 2 


V2 


-2 + c 2 


o 
z 


2 


-2 (256 + 304 c - 112 c 2 - 14 c 3 + c 4 ) 

(-2 + c) 2 (-8 + c 2 ) 


o 
z 


4 4- r- 2 


2 


2/2 2 


-3 (-104 + c 2 ) 

-8 + c 2 


2 V2 


-8 + c 2 


2 


2/i 2/3 


240 

-8 + c 2 


2 V2 


-8 + c 2 


3 


2/i 4 


1U + c 

6 (-2 + c) 


i 

2 




3 


2/i 2 2/2 


18 + c 

-2 + c 


1 


-1 + c 


3 


2/2 2 


16 

-2 + c 


2 


-2 + c 


3 


2/i 2/3 


22 + c 
-2 + c 


2 


-2 + c 
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V. THE CONNECTED PARTS 

The generating function of the connected parts of the 
average values of the total field reads 



with 



]n(R n (k))=a c nl k 2 +a c n2 k 4 + 



e (-i^«-^„(i)-i)!n^ 

i:Xift)=! m m ' 



(38) 



(39) 



We repeat that we are working exclusively in the HT 
phase and that we do not need to subtract powers of 
the magnetization. Using Eq. I|14l) and the construction 
discussed in the previous section, we can then calculate 
a n i(y»)- I n addition, we have 



2l\c 



(40) 



with the connected part of the average values () c defined 
in the usual way. For instance, 

<2 = an,2 - {l/2)a 2 n l 



2 n sites 



with 



{( 



E ^) 4 ) c = 

2 71 sites 

(( E ^') 4 > - 3<( 

2 ri s-i£es 



(41) 



E ^) 2 ) 2 

2" sites 



Wc define the finite volume susceptibility and their 
analog for the higher order derivatives of the free energy 
(zero momentum renormalized couplings) 



Xn 



(q) — ((Sg^sites fix)'') 



2 n 



(42) 



We restrict our considerations to the set of initial values 

(21) 

such that the infinite volume limit of \ n exists and 
is finite for every positive I. This means that we are 
not at another critical point or more generally not on a 
critical hypersurface at the boundary of the HT phase. 
We emphasize that the existence of the infinite volume 
limit requires < c < 2. For c > 2, the energy of a 
constant field configuration scales faster than the number 
of sites and the model has no interest from a statistical 
mechanics point of view. 

In the following, we assume that the initial values do,; 
are such that, 



lim x ( n ] = X 



(q) 



(43) 



is finite. From Eq. 1)40(1 . it is then clear that for n large 
enough, we have the leading scaling 



< ; « (2(-yr = Aft 



(44) 



It it thus tempting to find a simple relationship between 
a c n i and y n ,i- Indeed, such relation can be found at lowest 
non-trivial order from Eq. Q32[l which implies that 



(45) 



This can be seen either by using the Mobius inversion 
formula 

VI = E (-l) Io(iK1 x (X (i) - 1)! 

i:X x (i)=Z 



U 3 y? 



i 



III E t± ^r\ —v 



i r:Ii (r)— m 



or more simply by noticing that 



co ,21 



E^^nSr 



(47) 



Similar formulas are used in multiparticle scattering the- 
ory 00 

Eq. (|45|l means that there are no nonlinear contribu- 
tions of order /3 to af. For instance, there are no yf 
or yiy2 terms in ag. This is expected because the non- 
linear terms of order {3 l scale faster than yi, (assuming 
< c < 2). We we say that a term "scale faster", we 
mean that it goes to zero at a slower rate when n be- 
comes large. In general, at each RG step, a term fl m y^ 
of order (3 is multiplied by 



2 Xo(i) ( Si > A 
4 



(l) 



The strict inequality comes from the fact that for the 
nonlinear terms 2o(i) > E It is thus clear that nonlinear 
terms of order (3 l would spoil the HT scaling of Eq. I)44|) 
and contradict the existence of a infinite volume limit. 

For higher order terms, the sign of the denominator 
Di i introduced in Eq. (|27|) tells us whether or not the 
term scales faster or slower than the linear term. With 
our sign convention, c > c cr jt.(Z,i), means Di.i < and 
the term spoils the HT scaling Eq. H44|) . Since the co- 
efficients are rational functions of c, they cannot vanish 
suddenly when c becomes larger than c cr jt.(?,i). Con- 
sequently if < c cr jt.(Z,i) < 2, the coefficient of the 
corresponding term is expected to vanish identically. 

We have checked that this argument is consistent with 
our previous explicit calculations. We have used Eqs. 

(CHEjl and the already calculated coefficients in Eq. 
(1221 to calculate 



a i = vi + E ti ^v l y 

i:Xi(i)>2 



vl 2 



(48) 



up to order 7. For all the 50 terms with < c cr ;t. < 2, 
we found that the corresponding ti i are identically zero. 
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VI. THE HT PHASE 

In the previous section, we have argued (and checked 
explicitly up to order 7) that terms that scale faster than 
the linear term for c cr n. < c < 2 have a zero coefficient. 
In this section, we discuss more carefully some aspects of 
the argument and explain that having such terms nonzero 
would result in serious inconsistency. 

First of all, the existence of a HT phase is well es- 
tablished. The existence of a infinite volume limit j2^| 
and the absence of spontaneous magnetization for suffi- 
ciently high temperature 8] can be shown rigorously for 
< c < 2 and a Ising measure. Bounds on the free en- 
ergy density [2^, can be established for < c < 2 and 
measures with a compact support. The argument should 
also apply to measures that can be well approximated by 
measures with a compact support (see Eq. Q and the 
argument 26] that for Landau- Ginzburg measures, the 
restriction to |</>| < 4> m ax leads to exponentially control- 
lable errors). 

It is thus reasonable to assume that there exists some 
neighborhood of the HT fixed point where the infinite 
volume limit of the susceptibility and higher order deriva- 
tives (see Eq. (143(1 ) exist. Terms scaling faster than the 
linear term seem to contradict the existence of these in- 
finite volume limits. However, we should exclude the 
possibility that several terms (scaling identically) cancel 
each others. The existence of universal ratio of ampli- 
tudes means that we cannot in general pick arbitrary ini- 
tial values for the scaling variables. However, such con- 
straints apply for large values of the HT scaling variables. 
On the other hand, for arbitrarily small values of the HT 
scaling variables, one should be able to make indepen- 
dent variations of each variable while staying in the HT 
phase. This prevents the fine-tuning required to obtain 
cancellations. The HT fixed point R* = 1 corresponds to 
a local measure W{4>) oc S(4>) for which the correlations 
are zero. It is intuitively clear that by taking measures 
narrowly peaked at zero, one can avoid long range corre- 
lations. This continuity argument can probably be made 
rigorous by using Banach spaces as in Refs. [E E2- We 
conclude that the coefficients ti i in Eq. I|48|) of the terms 
with < c cr it.(l, i) < 2 must vanish identically. 



VII. THE ABSENCE OF POLES FOR < c < 2 

We are now in position to show that the small de- 
nominator problem can be evaded for any c such that 
< c < 2 and that the solution of the RG flows problem 
suggested in Eq. (|3*7}l can be constructed safely order by 
order. In section we have constructed the af in terms 
of the previously calculated a;. However we could have 
proceeded directly, writing the af t+1 ; in terms of the a c n f 

<+U = M k t a c nM + "?X,fc<« + ■■■ (49) 

k+q>l 



The coefficients v l q and the higher order ones can be 
obtained by using the expansion of Eq. (|38|l in the log- 
arithm of Eq. © and expanding order by order in a°. 
The series does not terminate. The linear transforma- 
tion is the same as before because af and a; only differ 
by nonlinear terms. Using 

<u = nK,r > (50) 

we obtain 

K+u = \(i)Kj + ]T wf q h c n ^ q + ... (5i) 

k+q>l 

We then introduce the expansion 

i:Xi(i)>/ m 

and obtain 



with given by a formula similar to Eq. 1(26(1 . except 
that it does not terminate. A detailed analysis shows that 
the two formulas have in common that the numerator 
depends only on coefficients of strictly lower orders in (3, 
and Eq. ((53(1 can be used order by order to construct the 
s^; for generic values of c. 

Since IZ^ 1 is upper triangular, we see from Eq. 1(50(1 
that hf is equal to af plus terms which go to zero faster. 
Consequently, for large n, the leading scaling is 

K,i « Afo • (54) 

Following reasonings used before, this implies that terms 
in the expansion Eq. 1(52(1 that scale faster than yi for 
any < c < 2 should have a vanishing coefficient. In 
other words: 

< Ccrtt.(M) < 2 =>■ sf ti = . 

Given the specific form of the sf ; given in Eq. 1(53(1 . the 
hf have no poles for < c < 2. The af being linear 
combinations of hf and the a/ being linear combinations 
of products of af, we conclude that the expansion of the 
ai in terms of the scaling variables have also no poles for 
< c < 2, in agreement with the conjecture stated in 
section Hvl 

Again we see that there exists a unique continuous def- 
inition of the scaling variables that can be used at par- 
ticular values of c where the denominator is exactly zero. 
From a practical point of view, the calculation at fixed 
c of the sf j is easier than the calculation of the sj i, be- 
cause no limit needs to be taken explicitly. The sf ; being 
rational function of c they cannot be zero everywhere ex- 
cept at isolated values. Consequently, we can set to zero 
the sf j having c cr u < 2 even at values of c where 
Dia = 0. 
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VIII. THE INITIAL VALUE PROBLEM 



IX. THE HT EXPANSION 



We now return to Eq. 1(3*71) . In order to complete 
our solution of the problem, namely expressing the a„ 
in terms of their initial values ao, we need to calculate 

y(a )- 

Before doing this, we want to show that the initial 
values yo have a very simple interpretation. We have 
learned in the previous sections that y n j is the only 
leading term of a c n l when n becomes large. If at a given 
< c < 2, a nonlinear terms scales exactly like y n ,i, then 
by increasing c slightly (but keeping c < 2), we can make 
this term dominant in contradiction with the existence 
of the infinite volume limit. Consequently, 



lim A ; n a c n t — lim A ; 



From Eq. (|4T"|) , we see that 



Un,l = yo, I 



(21) 



(55) 



(56) 



This means that the infinite volume limit of the suscepti- 
bility and of the the higher derivatives of the free energy 
density completely determine the RG flows in the HT 
phase. This also means that the calculation of yo,; given 
ao.; is nontrivial far away from the HT fixed point. How- 
ever, we can take advantage of the fact that 



\i)Vn,i = yo,, 



(57) 



to estimate yo,z using expansions valid at intermediate 
values of n. 

We now discuss the inversion question. We need to 
determine the coefficients ry of the expansion 



A simple application of the method presented here is 
the calculation of the high-temperature expansion at fi- 
nite volume. As a simple example, we consider the first 
order correction to the susceptibility for a Ising measure 
(Ro(k) = cos(y/j3k)). Using the results found in the pre- 
vious sections, we obtain 



X 



(2) 



-2 „ j^\n 







6c c 

+- — (7) yo,2 



4- c v 4 
1 +(3b 1 . n + 0((3 2 ) 



Using ao, 



if- and ao,2 



we obtain 



2/o,i = - 

2/0,2 = - 

and consequently 

h,n = 



f3 2 c 



I3 2 c{2 + c) 



2 4(4 -c) 4(4-c)(2-c) 
12 ' 



2c —<-Y 

(4-c)(2-c) 2-c2 ] 



4-c 4' 



(61) 



(62) 



This is in agreement with results obtained [19| using 
graphical methods. 



Vl 



(58) 



X. RG INVARIANTS 



This can be done by replacing the yi appearing in the 
expansion of the hi in Eq. by Eq. . This yields 
an Eq. of the form 

n,i + ay + X Li = , 

with Xn linear in the s and multilinear in r of strictly 
lower order. One can then construct the ry order by 
order without ever creating a pole in the range < c < 2. 
At lowest non-trivial order, we have 



h l= E (-l) Io(iK1 (Io(i)-l)!n^V+°(/ 3 ' +1 ) 



i:Z 1 (l)=i 



(59) 



A more detailed analysis shows that for higher orders 

c 



ry - (^-) Il(i) -^,i(c) , 
c — 4 



(60) 



with Yi.i(c) having poles only for 2 < c < 4. The values 
of Yi.i{c) up to order 4 are given in Table ITT1 



In Hamiltonian mechanics, integrable systems with q 
degrees of freedom have q constants of motions and q 
periodic variables with independent periods depending 
on the constants of motion. In the present case, time 
is discrete and exponential decays replace the quasiperi- 
odic behavior. For a truncation of dimension l m axi it is 



nevertheless possible to construct l Tl 
motion: 



Gi 



-(21) 



Un,l 



(-2» n ,i)('-D(o/2)+i 



1 constants of 



(63) 



These quantities are n-independent and we call them RG 
invariants. We can calculate them at n = 0. Using Eq. 
(|56|l . we obtain 



Gi - (-1) 



i+i 



( x (2))(/-l)(D/2)+Z 



(64) 



We can also calculate them at large enough values of n 
where the HT expansion works well. The minus sign has 
been introduced in order to have Gi > for D = 3. 
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TABLE II: Values of Yi,i(c), and T t ,i(c) defined in the text. 



h 1 3 
h 1 h 2 



h! 2 h 2 



hih 3 



-m 



160-224c+44c 2 +4c 3 +c 4 
(-2+c) 2 (-8+c 2 ) 

3(-40+c 2 ) 
-8+c 2 



-229376+638976c-364544c 2 -160768c 3 + 109056c 4 +11648e 5 + 1664c 6 -5952c 7 +392e 8 + 148e 9 +6e 10 +e 11 
(-2+c) 3 (-8+c 2 )(-32+c 3 )(-16+c 3 ) 

-6(24576~53248c- 19456c 2 +10496c 3 +3968c 4 +800c s -576c 6 -40c 7 + 2c 8 +c 9 ) 
(-2+c)(-8+c 2 )(-32+c 3 )(-16+c 3 ) 

-1536 

-32+c 3 

-15 



/life 

h l 2 h 2 

h 2 2 
h\hz 



6+c 
-2+c 

_ ( 14+c N 
1,-2+cj 

- (608- 1504c+356c 2 +36c 3 +3c 4 ) 
2(-2+c) 2 (-8+c 2 ) 

-2(-704+1216c-316c 2 -20c 3 +c 4 ) 
(-2+c) 2 (-8+c 2 ) 

3(-104+c 2 ) 
-8+c 2 

-240 



-8+c 2 



ftl 4 

ftl 2 /l2 
h 2 2 



( 10 + c N 
1,-2 + cj 



2(18+c) 
-2+c 



-2+c 



f 22+c N 
1,-2+cj 



We now concentrate on the unstable direction of Wil- 
son's fixed point. We set the relevant scaling variable 
associated with this fixed point to a value u which be- 
comes our coordinate along the unstable direction and 
we set all the irrelevant ones to zero. We call Gi(u) the 



corresponding value of the ratio. Given that u is a scaling 
variable and that Gi is RG invariant, we have 

Gi(X w u) = Gi(u) , (65) 

with X w the eigenvalue corresponding to the unstable di- 
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rection of Wilson's fixed point. Consequently, we have 
the Fourier expansion: 

G l (u) = J2ArU tr " , (66) 

r 

with 



If the oscillatory terms are very small, as noticed in Refs. 
[HiL I27I l28| , we have the approximate universal ratios 

G,( M )~A ; , . (68) 

These constants can be calculated in an intermediate re- 
gion where the expansions in both scaling variables are 
valid 0. 

XI. CONCLUSIONS 

We have shown that the scaling variables correspond- 
ing to the HT fixed point of Dyson's hierarchical model 
can be constructed order by order without small denom- 
inator problems. The ambiguity noticed before ^3 f° r 
calculations at fixed values of c can be raised by requiring 
the continuity in c. Practical calculations at finite c are 
most easily done by following the explicit construction 
sketched in section fVlII for the connected part where no 



complicated limit is required. The remaining poles for 
2 < c < 4 reflect the degeneracy of the linear spectrum 
at c = 4 or the fact that some finite size corrections be- 
come leading effects for some value of 2 < c < 4 (where 
the infinite volume limit does not exist). 

We have solved the linear problem in compact form 
but at this point no compact form is available for the 
nonlinear problem. Even though we have "constants of 
motion" (the RG invariants), we do not have simple ex- 
pressions for them (as the integrals of integrable models). 
The question of the convergence of the series remains to 
be addressed and should result in the construction of the 
boundary of the HT phase. Finally, it would be desir- 
able to extend the method to models with 77 7^ 0. One 
way to achieve this goal would be to develop methods to 
systematically improve the hierarchical approximation. 
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